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tory  at  Rock  Island  under  an  In-House  Laboratory  Independent  Research  Project.  The 
purpose  of  the  investigation  was  to  analytically  study  the  thermal  radiation  heat 
transfer  inside  gun  tubes.  A  simplified  physical  model  was  used  in  the  analysis. 

With  this  model,  the  steady  turbulent  flow  of  an  optically  thick  gray  gas  in  the  gun 
tube  was  considered.  Both  developing  and  fully  developed  turbulent  flow  were  con¬ 
sidered.  A  mathematical  model  was  derived  and  programmed  for  numerical  evaluation. 
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heat  transfer.  The  absorption  coefficient  of  the  propellant  gas  is  very  high  because 
of  high  pressure  and  density.  Thus,  while  a  significant  amount  of  radiant  energy 
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ABSTRACT 


This  investigation  was  conducted  by  the  Research  Directorate  of  the 
Weapons  Laboratory  at  Rock  Island  under  an  In-House  Laboratory  Indepen¬ 
dent  Research  Project.  The  purpose  of  the  investigation  was  to  analyti¬ 
cally  study  the  thermal  radiation  heat  transfer  inside  gun  tubes.  A 
simplified  physical  model  was  used  in  the  analysis.  With  this  model, 
the  steady  turbulent  flow  of  an  optically  thick  gray  gas  in  the  gun 
tube  was  considered.  Both  developing  and  fully  developed  turbulent 
flow  were  considered.  A  mathematical  model  was  derived  and  programmed 
for  numerical  evaluation.  The  numerical  results  were  plotted  in  dimen¬ 
sionless  form  so  that  they  could  be  used  to  estimate  the  radiation  heat 
transfer  in  any  gun  tube  for  which  the  propellant  gas  flow  conditions 
are  known.  The  contribution  of  radiative  heat  transfer  in  the  XM140 
gun  tube  was  examined.  Radiation  heat  transfer  was  significant  near  the 
breech.  At  locations  removed  from  the  breech,  turbulent  convection  was 
the  dominant  mode  of  heat  transfer.  The  absorption  coefficient  of  the 
propellant  gas  is  very  high  because  of  high  pressure  and  density.  Thus, 
while  a  significant  amount  of  radiant  energy  is  emitted  by  the  propel¬ 
lant  gas,  most  of  this  energy  is  absorbed  before  it  strikes  the  bore 
surface.  Radiation  is  more  important  in  weapons  operating  at  lower 
pressures  and  having  lower  projectile  velocities. 
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nomenclature 


Symbols 

c  =  Specific  heat 
eb  =  Black  body  emission 
f  ~  Friction  factor 
9C  ~  Gravitational  constant 
N  -  Heat  transfer  coefficient 
K  =  Thermal  conductivity 
k  *  Absorption  coefficient 

kR  -  Rosseland  mean  absorption  coefficient 

i  -  Length 

• 

m  =  Mass  flow  rate 

Nu  =  Nusselt  number 

P  '■=  Pressure 

Pr  -  Prandtl  number 

0  *  Heat  flux 

r  =  Radial  coordinate 

rQ  -  Outside  radius 

ro*  -  Dimensionless  outside  radius 

Re  s  Reynolds  number 

T  '  Temperature 

t  =  T ime 

u  =  Velocity 

u*  -  Dimensionless  velocity 
v  =  Velocity 
Vp  =  Projectile  velocity 
x  =  Spat’ al  coordinate 
y  ~  Spatial  coordinate 
y*  =  Dimensionless  coordinate 
3  “  Thermal  expansion  coefficient 

5  =  Boundary  layer  thickness 

e  -  Eddy  diffusivity,  emissivity 

3  ~  Dimensionless  temperature 


iv 


X  =  Wavelength 

m  =  Viscosity 

v  =  Viscosity 

5  =  Dimensionless  parameter 

p  =  Density 

o  =  Radiation  parameter 

t  =  Optical  length 

*  =  Rayleigh  dissipation  function 

Subscripts 

d  =  Based  on  diameter 

9  =  Gas  value 

R  =  Radiation  value 

W  =  Wall  value 

X  =  Based  on  distance 

X  -  Wavelength 

o  -  Initial  value 

00  =  Free  stream  value 
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INTRODUCTION  '  ■ 

t 

The  purpose  of  tlie  present  investigation  is  to  analytically  inves¬ 
tigate  the  contribution  of  thermal  radiation  to  gun  tube  heating  when  a 
hot  propellant  gas  flows  through  the  gun  tube.  Thermal  radiation  has 
been  disregarded  in  previous  gun  tube  heat  transfer  analyses.  This  dis¬ 
regard  was  justified  on  the  basis  of  the  use  of  the  following  equation:1 

»  *  i  ,  ■ 

< 

'  q  =  0£g€w  (Tg-  -  V)  .  '  '  0) 

f 

The  maximum  radiative  heat  f.ux  was  calculated  by  the  assignment  of 
values  of  unity  to  eq  and  t  .'  The  calculated  radiative  flux  was  shown 
to  be  much  less  thanythat  of  experimental  values.  This  fact  was  used 
as  a  justification  for  disregarding  thermal . radiation.  Unfortunately, 
Equation  1  does  not  apply  to  the  physical  situation  encountered  in  gun 
tubes.  Equation  1  is  valid  for  only  an  isothermal  homogeneous  gas  sur¬ 
rounded  by  ah  isothermal  surface  and  in  the  absence  of  simultaneous 
conduction  or  convection.  The  propellant  gas  may  be  characterized  as 
optically  thick.'  Radiation,  in  an  optically  thick  medium  may  be  described 
as  a  diffusion  process  having  a  "radiation  conductivity"  given  by2 

.  '  k  =  IMi  1  '  ‘  (2) 

R  3k  1 

With’ the  use  of  Equation  2,  under  certain  circumstances,  the  radiation 
conductivity  in  propellant  gases1  is  much  greater  than  the  thermal  con¬ 
ductivity.  Thus,  the  conclusion  might  be  drawn,  erroneously,  that  radi¬ 
ation  i.s  the  dominant  mode  of  heat  transfer  and  that  convection  may  be 
disregarded.  Both  of  the  approaches  described  are  oversimplified 
since  the  essential  character  of  energy  transfer  in  a  medium  which  ab¬ 
sorbs  and  emits  radiation  Ts  disregarded.  The  total  heat  transfer 
mechanism,  vyhich. includes  the  interaction  of  conduction,  convection,  and 
■radiation,  must  be  considered  in  any  rigorous'  analysis.  The  study  of 
the  interaction  o*  radiation  with  conduction  and  convection  in  ab¬ 
sorbing  and  emitting  media  is  a  relatively  new  field.  The  greater  part 
of  this: work  has  been  conducted  since  1960.  V iskanta3’4  studied  the 
interaction  of  laminar  convection,  and  radiatibn  in  channels  and  tubes. 
Einstein5  considered  radiant  heat  transfer  to  absorbing  gases  inclosed 
in  a  circular  pipe  with  conduction,  gas  flow,  and  internal  heat  genera¬ 
tion.  Viskanta6  has  also  compiled  an  extensive  review  on  the  subject. 
More  recently,  Pearce  and  Emery7  studied  the  combined  heat  transfer  by 
radiation  and  laminar  convection  to  a  radiatibn  absorbing  fluid  in  the 
entry  region  of  a  pipe.  Substantially  less  work  has  been  done  on  the 
inter„ -tion  of  radiation  and  turbulent  'flows.  Nichols8  investigated 
radiation  in  the  turbulent  flow  of  steam  in  the  entrance  region  of  an 
annulus.  Landram  et  a!-  analytically  investigated  heat  transfer  in 
turbulent  pipe  few  with  optically  thin  radiation. 
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The  propellant  gas  flow  in  a  gun  tube  is  unsteady  and  turbulent. 

The  analysis  of  unsteady  turbulent  flow  is  a  new  field  and  not  well  ad¬ 
vanced.  Heat  transfer  analyses,  for  this  type  of  flow  usually  involve 
solutions  of  the  transient  -fluid  dynamics  equations  and,  then,  use  of 
some  analogy  to  determine  the  heat  transfer  coefficient.10  No  refer¬ 
ences  were  found  in  which  an  attempt  was  made  to  compute  the  unsteady 
turbulent  temperature  field.  Knowledge  of  the  temperature  field  is 
required  to  determine  the  interaction  of  convection  with  radiation 
since  the  radiant  heat  flux  is  a  function  of  the  local  temperature. 

Thus  the  calculation  of  radiant  heat  flux  in  unsteady  turbulent  flow 
such  as  is  found  in  gun  tubes  is  beyond  the  present  state  of  the  art. 

The  approach  in  the  present  investigation  is  to  consider  a  constant 
property,  gray  gas  (one  whose  radiation  properties  are  independent  of 
wavelength)  for  steady  turbulent  flow  in  a  tube.  While  this  is  an  ex¬ 
treme  simplification,  the  present  approach  should  provide  a  substan¬ 
tially  better  method  for  predicting  the  radiative  flux  in  gun  tubes 
than  any  other  previous  method. 

A  discussion  of  radiation  in  participating  media  is  given  in  the 
following  section.  The  analytical  section  deals  with  the  analyses  of  a 
constant  property,  gray,  optically  thick  gas  for  steady  undeveloped 
turbulent  flow  and  for  steady  fully  developed  flow  in  a  tube.  The 
analytical  results  are  presented  and  discussed, 

RADIATION  IN  PARTICIPATING  MEDIA 

The  expression  for  the  conservation  of  energy  in  a  moving  fluid  is 


r  Of  - 
pCp  Dt  ' 


div  q 


61  K  *  -♦ 


Equation  3  is  equally  applicable  to  both  radiation  participating  and 
nonparticipating  fluids.  The  difference  is  in  the  expression  for  the 
heat  flux  vector,  q.  When  both  conduction  and  radiation  are  present, 
the  total  heat  flux  vector  may  be  expressed  as 

q  *  -  K  grad  T  t  qR  ^ 

where  -  K  grad  T  expresses  thermal  conduction  and  qR  represents  thermal 
radiation.  The  radiation  flux  vector  mus1,  be  formulated  as  a  function 
of  temperature  to  solve  Equation  4.  The  heat  flux  vector  can  be  an 
extremely  complicated  mathematical  function.  This  vector  is  dependent 
upon  the  flow  field,  the  gas  properties,  and  the  geometry  of  the  sur¬ 
rounding  inclosure  To  reduce  the  heat  flux  vector  to  a  mathematically 
tractable  form  usually  requires  a  number  of  simplification  assumptions. 
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One  of  the  most  important  parameters  by  which  radiation  is  described 
in  a  participating  medium  is  the  optical  thickness,  t.. 

A 


T 


A 


(5) 


The  optical  thickness  represents  the  ratio  of  a  characteristic  dimension 
of  the  system  to  the  photon  mean  free  path  length.  The  value  of  the 
optical  distance  characterizes  three  separate  radiation  regimes.  When 
the  value  of  t.  is  much  less  than  one,  the  radiation  is  called  optically 
thin.  The  photon  mean  free  path  is  much  greater  than  the  system  dimen¬ 
sions.  Any  photon  emitted  by  a  gas  volume  will  travel  directly  to  the 
bounding  surfaces  without  interacting  with  any  other  gas  volume.  Thus, 
the  fluid  absorbs  only  radiation  emitted  by  the  boundary  surfaces.  This 
is  also  referred  to  as  negligible  self-absorption.  When  •.  is  approxi¬ 
mately  equal  to  one,  radiation  exchange  occurs  between  al'rthe  fluid 
elements,  and  between  all  fluid  elements  and  the  boundaries.  This  re¬ 
sults  in  an  integral  expression  for  the  radiation  flux  vector. 

The  third  regime  involves  situations  where  t  is  much  greater  than 
one.  The  photon  mean  free  path  length  is  much  less  than  the  system  di¬ 
mensions  and  the  radiation  may  be  characterized  as  a  photon  continuum. 
Radiation  transfer  within  the  medium  becomes  a  diffusion  process.  As 
shown  by  Sparrow  and  Cess,2  the  heat  flux  vector  for  this  regime  is 
given  by 


q 


RA 


3^  grad  e 


bA 


16) 


Equation  6  is  valid  within  one  photon  mean  free  path  length  from  the 
bounding  surface.  Beyond  a  distance  of  one  mean  free  path  length,  no 
photons  emitted  b>  the  surface  will  be  present.  Thus,  radiation  from  the 
boundary  will  not  be  detected  by  the  medium,  and  the  heat  transfer  within 
the  medium  will  be  independent  of  the  surface  emittance  Some  question 
exists  as  to  whether  Equation  6  applies  at  the  boundary  itself.  However, 
Sparrow  end  Cess  argue  that  if  Equation  6  were  not  valid  at  the  boundary, 
a  discontinuity  in  radiation  flux  would  occur,  which  is  physically  im¬ 
possible.  They  conclude  that  Equation  6  is  valid  at  the  boundaries  and 
throughout  the  medium  when  optically  thick  conditions  prevail. 


Equation  6  gives  the  heat  flux  vector  as  a  function  of  wavelength. 
The  assumption  of  a  gray  gas  eliminates  the  dependence  on  wavelength. 
The  total  radiation  heat  flux  is  determined  by  integration  over  all 
wavelengths 

JL  debx  Hi  /7\ 

3kx  dy  dx  (7) 


^  R  = 


3 


Integration  of  Equation  7  results  in 


.4  deb  16oT-  cP 
qR  "  3kR  dy  "  3kR  dy 


where 


k 


(8) 


(9j 


The  quantity  kR  is  called  the  Rosseland  mean  absorption  coefficient, 
and  the  quantity  16'oF-73kR  is  the  radiation  conductivity 

The  absorption  coefficient,  ,  must  be  known  to  evaluate  Equation 
9.  The  function  k  has  been  experimentally  determined  tor  a  number  of 
pure  gases  at  room" temperature  and  at  atmospheric  pressure  The  absorp¬ 
tion  or  emission  of  radiation  by  a  gas  occurs  at  rued  frequencies  over 
narrow  frequency  bands.  The  emission  and  absorption  at  these  frequen¬ 
cies  are  due  to  transitions  between  energy  levels  of  the  atoms  and 
molecules  that  constitute  the  gas  The  frequency  bands  may  be  collision 
broadened.  This  phenomenon  occurs  when  the  pressure  and  temperature 
become  high,  and  the  frequency  of  molecular  and  atomic  collisions  is 
increased.  In  mixtures  of  gases,  an  overlapping  or  absorption  bands 
may  exist.  While  k  may,  in  theo-y,  be  computed  nom  quantum  mechanics, 
experimental  data  are  usually  ne^essa  y  unfortunately,  the  locating 
of  any  experimental  data  tor  the  absorpt  on  v,oefnc»ent  us  a  function 
of  temperature  and  pressure;  and  tor  the  currently  used  propellant  gases 
was  impossible.  The  prima-y  radiation  participating  constituents  of 
propellant  gases  are--  carbon  monoxide,  54  percent  by  weight;  ca-bon 
dioxide,  17  percent  by  weight;  and  water  vapor,  15  percent  by  weight. 
Tien1*  has  used  experimental  data  to  compute  the  Rosseland  mean 
absorption  coefficient  for  each  of  these  gases  as  functions  of  tempera¬ 
ture  and  pressure.  To  simply  add  the  absorption  coefficient  ot  each 
constituent  gas  to  find  the  total  absorption  coefficient  is  impossible 
since  some  of  the  absorption  bands  overlap  However,  to  obtain  the 
minimum  value  of  absorption  coefficient  is  possible  since  it  is  at 
least  as  large  as  the  largest  absorption  coefficient  oi  any  of  the 
constituent  gases-  This  minimum  mean  aosorpnon  coefr'c-ent  can  os 
combined  with  some  characteristic  length  to  determine  a  characteristic 
optical  dimension  A  parametric  study  can  then  be  performed  Dy  variation 
of  the  optical  dimension  The  results  will  show  the  range  of  optical 
distances  for  which  radiation  in  gun  tubes  is  important.  Thus,  rational 
judgments  can  be  made  about  the  contribution  of  radiation  to  Overall 
gun  tube  heating 
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The  approximate  minimum  value  of  kR  may  be  determined  from  Figures 
22,  23,  and  24  of  Tiens  article,  and  is 


,  „  60 

kR  .  “  atm  ft 

nnn 


(10) 


For  a  gun  tube  pressure  of  30,000  psi  and  a  bore  diameter  of  0.3  inch, 
the  optical  diameter  has  a  value  of  3000.  This  value  is  much  greater 
than  one,  and  the  radiation  is  in  the  optically  thick  regime. 

THEORETICAL  ANALYSIS 


The  analysis  of  radiation-convection  interaction  for  optically 
thick  turbulent  flow  in  a  tube  is  outlined  in  this  section.  The  assump¬ 
tion  is  that  fluid  properties  are  constant;  therefore,  the  fluid  flow 
is  independent  of  the  temperature  distribution.  Additional  assumptions 
are  that  the  flow  is  steady  and  that  viscous  dissipation  is  negligible. 
The  flow  of  both  undeveloped  and  fully  developed  turbulence  is  consid¬ 
ered.  The  various  flow  regimes  for  turbulent  flow  in  a  tube  are  illus¬ 
trated  in  Figure  1. 

Undeveloped  Turbulent  Flow 

When  the  thickness  of  the  hydrodynamic  boundary  layer  is  much  less 
than  the  tube  radius,  the  flow  is  not  fully  developed  and  is  closely 
approximated  by  the  flow  of  a  semi -infinite  fluid  over  a  flat  plate. 

The  thickness  of  the  thermal  boundary  layer  is  nearly  the  same  as  the 
thickness  of  the  hydrodynamic  boundary  layer  if  the  Prandtl  number, 

Pr  =  v/a,  is  close  to  unity.  This  condition  is  true  for  most  gases. 

The  equality  of  the  boundary  layer  thicknesses  will  be  assumed  in  this 
analysis.  Equation  3,  subject  to  all  the  restrictions  mentioned,  may 
be  written  as 


3T  _  3  p  /  k  16ol  3 
3y  ~  ay  u  pc  T  3kRpc 


)  3 

'  ay  J 


01) 


where  e  is  the  turbulent  eddy  diffusivity  for  heat  transfer.  The  case 
to  be  considered  is  that  in  which  the  tube  walls  are  at  a  constant  tem¬ 
perature,  To,  and  the  inviscid  core  is  at  a  constant  temperature,  Tc^ 

A  three-layer  model  will  be  used  to  describe  the  turbulent  flow.  The 
three  layers  are  the  laminar  sublayer,  buffer  layer,  and  turbulent 
layer.13  In  the  two  sublayers  near  the  wall,  the  radial  velocity,  v, 
is  approximately  equal  to  zero.  Since  the  tube  surface  temperature  is 
constant,  3T/3x  is  also  approximately  equal  to  zero  in  the  sublayers. 
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The  approximate  minimum  value  of  kR  may  be  determined  from  figures 
22,  23,  and  24  of  Tiens  article,  and  is 

=  ItmTt  (10) 


For  a  gun  tube  pressure  of  30,000  psi  and  a  bore  diameter  of  0.3  inch, 
the  optical  diameter  has  a  value  of  3000.  This  value  is  much  greater 
than  one,  and  the  radiation  is  in  the  optically  thick  regime. 

THEORETICAL  ANALYSIS 

The  analysis  of  radiation-convection  interaction  for  optically 
thick  turbulent  flow  in  a  tube  is  outlined  in  this  section.  The  assump¬ 
tion  is  that  fluid  properties  are  constant;  therefore,  the  fluid  flow 
is  independent  of  the  temperature  distribution.  Additional  assumptions 
are  that  the  flow  is  steady  and  that  viscous  dissipation  is  negligible. 
The  flow  of  both  undeveloped  and  fully  developed  turbulence  is  consid¬ 
ered.  The  various  flow  regimes  for  turbulent  flow  in  a  tube  are  illus¬ 
trated  in  Figure  1. 

Undeveloped  Turbulent  Flow 

When  the  thickness  of  the  hydrodynamic  boundary  layer  is  much  less 
than  the  tube  radius,  the  flow  is  not  fully  developed  and  is  closely 
approximated  by  the  flow  of  a  semi-infinite  fluid  over  a  flat  plate. 

The  thickness  of  the  thermal  boundary  layer  is  nearly  the  same  as  the 
thickness  of  the  hydrodynamic  boundary  layer  if  the  Prandtl  number, 

Pr  =  v/a,is  close  to  unity.  This  condition  is  true  for  most  gases. 

The  equality  of  the  boundary  layer  thicknesses  will  be  assumed  in  this 
analysis.  Equation  3,  subject  to  all  the  restrictions  mentioned,  may 
be  written  as 


16qT3 

3k^pc 


+  i-  )  |I  ] 
'  ay  J 


where  e  is  the  turbulent  eddy  diffusivity  for  heat  transfer.  The  case 
to  be  considered  is  that  in  which  the  tube  walls  are  at  a  constant  tem¬ 
perature,  To,  and  the  inviscid  core  is  at  a  constant  temperature,  T*. 

A  three-layer  model  will  be  used  to  describe  the  turbulent  flow.  The 
three  layers  are  the  laminar  sublayer,  buffer  layer,  and  turbulent 
layer.13  In  the  two  sublayers  near  the  wall,  the  radial  velocity,  v, 
is  approximately  equal  to  zero.  Since  the  tube  surface  temperature  is 
constant,  3T/3x  is  also  approximately  equal  to  zero  in  the  sublayers. 
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I 


Thus,  for  the  sublayers,  Equation  11  reduces,  to 


^  K  x  +  16211 ;  £ )  f  .  0 

,dy  1 ■  pc  ’  3k^pc  ,  •  '  dy 


Equation  12  may  be  integrated  once  to  give 

■{  K  ,  ,  16aT3  >  .  ,  dT 

'  (  W  +  3iyc  +  £,  )  dy'  -  c  • 


When  y  =  0,  e  =  0  arid  , '  > 

!  ‘  I 

{  JL  4.  16afo3  v  dT  _  n. 

1  /  PC:  3kRpc  '  dy  ~  "q' 


where  q  >  0. indicates  that  heat  is  being  added  to  the  fluid.  Equation 
14,  is  used  to  find  'the  constant  of  integration  and  Equation  13  may  be 
rearranged  to  give  1  ! 

1  • 

'  dT  _  -  "V  ,, _ ,  ’  ’ 

dy  "  77  K1  .  I60T3  .  •.  ,  ■  (15) 

.  •  pG  (  7c  r7k^  ~  e  }  *.  1 


,  ,  A  differed  procedure  is  used  to  determine  the'  temperature  distri¬ 
bution  in  the  turbulent  layer,  Reynolds analogy  is  used  in  the  turbu¬ 
lent  layer.  The  relation  between  the  velocity  and  temperature  profiles 


du  = 


Equation  16  is  valid,  only  if  the  total  di f fus i vi ties  for  heat  and  momen¬ 
tum  are  equal  and  if  the  turbulent  diffusi vities  are  much  greater  than 
the  molecular,diffusivities.  These  conditions  are  usually  satisfied  in 
the  turbulent  flow  of  gases  when  no  radiation  is  present.  Some  care 
must  be  exercised  in  the  use  of  Equation  16  for  the  present  investiga¬ 
tion.  , The  total  thermal  conductivity  is  the  sum  of  the  molecular  con¬ 
ductivity  and  the'  radiation  conductivity.  If  the  radiation  conductivity 
is  much  higher  than  the  molecular  conductivity,  the  Prandtl  number  cal¬ 
culated  from  the  total  conductivity  may  become  quite  small.  Molecular 
diffusion  may  be  of  the  same  order  of  magnitude  as  turbulent  diffusion 
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in  this  case.  Equation  16  will  be  valid  only  if  the  quantity  s/v  is 
much  larger  than  the  quantity  (  1  +  16cT3/3kRK)/Pr  throughout  the  tur¬ 
bulent  layer.  This  requirement  may  be  difficult  to  satisfy  near  the  edge 
of  the  buffer  layer. 

Experimental  data  for  the  velocity  profiles  in  steady  turbulent  flow 
have  been  correlated  so  that  a  dimensionless  velocity,  u*,  is  a  known 
function  of  a  dimensionless  coordinate,  y*.  This  function  is  referred 
to  as  the  universal  velocity  profile  for  turbulent  flow.  To  utilize 
these  data,  the  governing  equations  must  be  expressed  in  terms  of  these 
dimensionless  variables,  which  are  defined  as 


/  Vo/p 


(17) 


and 


y*  =  y 


/o  T  . 

'C  o/p 


(18) 


Also,  the  expressing  of  temperature  in  a  dimensionless  form  is  desirable 
so  that  all  significant  parametric  groups  affecting  the  temperature  dis¬ 
tribution  can  be  isolated.  The  following  dimensionless  temperature  is 
defined: 


e 


T  ( 


*  )  Pr  Re 
o 


(19) 


The  introduction  of  the  dimensionless  variables  into  the  differential 
equation  for  the  sublayers  gives 


de_ 

dy* 


-  1 


(%  (pf 


4  0' 
sx 


Pr4Re* 

A 


£1 


where 

Sx 


(20) 


(21) 


8 


The  friction  factor  for  steady  turbulent  flow  over  a  flat  plate  is13 


fv  =  0.059  Re  ~0,2 

A  X 


(22) 


Thus,  the  governing  differential  equation  for  the  sublayers  becomes 


de 

dy* 


[ 


-  5.822  ReY 
Pr  +  Pr4Re3 

A 


0.1 


(23) 


The  dimensionless  form  of  the  differential  equation  for  the  turbulent 
layer  is 


de 

du* 


-  5.822  Re 

X 


o.i 


(24) 


Several  different  functional  forms  for  the  universal  velocity  pro¬ 
file  and  eddy  diffusivity  have  been  proposed.  Kays13  recommends  the 
following  relations: 


Laminar  Sublayer  0  <  y*  <  5 


u* 


(25) 


Buffer  Layev*  5  «•  y*  <  30 


u* 


-  3.05  +  5  In  y*  ,  J  - 


(26) 


Equation  24  may  he  integrated  directly  to  give 


e  =  5.822  u*  Re"0,1  +  c  (27) 

A 


The  constant  of  integration  is  determined  by  evaluation  of  Equation  27 
at  y*  =  30  where,  from  Equation  26,  u*  =  14.  The  dimensionless 
temperature  distribution  in  the  turbulent  layer  is  given  by 


0(y*)  =  0(30)  -  5.822Rev0el  (u*  -  14)  (28) 

A 

When  Equation  28  is  evaluated  in  the  free  stream,  y*  =  «®  and  U*~  = 
l/^fx/2,  the  following  result  is  obtained: 

0(»)  -  0(30)  =  -  5.822Rev0,1  (5.822Rev0,1  -14)  (29) 

X  A 

The  set  of  equations  from  which  the  total  temperature  difference 
across  the  turbulent  boundary  layer  can  be  obtained  are 

Laminar  Sublayer,  0  <  y*  <  5 


de 

3y* 


5.822Rex°'^  [ - ^3-] 

1  +  - X - 

(PrRex)3 


Buffer  Layer;  5  <  y*  <  30 


de 

dy* 


5.822Rev0,1 

A 

€x03  ,  y* 

+  5 

A 


-I  3 


Turbulent  Layer,  30  <  y*  <  ~ 

e(»)  -  0(30)  =  -  5.822Re  °* 1  (5.822Rev0,1  -14) 

A  A 


(30) 


The  heat  transfer  coefficient  is  defined  as 


h  = 


_9°_ 


(VTo) 

and  the  Nusselt  number  is  defined  as 
Nu„ 


hx 

K 


gox 

K(T  -To) 


(31) 


(32) 
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When  the  dimensionless  variables  are  introduced  into  Equation  32,  the 
following  expression  for  the  Nusselt  number  is  obtained: 


Pr  Rex 

Nux  ”  eR  -  e(o) 


(33) 


The  set  of  governing  equations  for  the  temperature  profile  may  be 
solved  numerically  by  the  fourth  order  Runge-Kutta  method.14  The 
Nusselt  number  is  then  evaluated  from  Equation  33.  The  parameter  5X 
may  be  varied  to  give  the  Nusselt  number  as  functions  of  the  Reynolds 
number  and  ?x.  These  results  will  apply  for  the  case  where  the  heat 
flux  at  the  tube  wall  is  specified.  A  different  procedure  of  solution 
is  required  when  the  wall  temperature  and  free  stream  temperature  are 
specified.  The  given  wall  temperature  is  specified  as  the  initial 
value.  The  desired  values  are  assigned  to  Rex,  Pr,  (ov/Ku<J,  and 
(ucoArv).  The  value  of  (qgX/K)  is  adjusted,  through  an  iterative  pro¬ 
cedure,  until  a  solution  with  the  correct  free  stream  temperature  is 
obtained.  The  results  give  the  Nusselt  number  as  a  function  of  the 
Reynolds  number,  the  Prandtl  number,  wall  temperature,  free  stream  tem¬ 
perature,  thermal  conductivity,  and  (u^/kpv).  The  results  of  a  numeri¬ 
cal  evaluation  are  given  in  the  following  section. 


Fully  Developed  Turbulent  Flow 

The  flow  in  a  tube  is  fully  developed,  in  the  hydrodynamic  sense, 
when  the  thickness  of  the  hydrodynamic  boundary  layer  is  equal  to  the 
tube  radius.  The  velocity  profile  in  the  tube  is  then  independent  of 
axial  location,  and  the  radial  velocity  component  is  zero.  However, 
this  does  not  imply  that  the  flow  is  fully  developed  in  the  thermal 
sense  since  the  fluid  temperature  will  still  be  a  function  of  the 
axial  location.  The  temperature  gradient  in  the  axial  direction  will 
be  a  constant  if  the  heat  flux  at  the  tube  wall  is  a  constant.  Equation 
3,  subject  to  all  previous  assumptions  and  the  additional  assumptions 
of  fully  developed  flow  and  constant  heat  flux,  is 


u  41  =  LA.  r  r  ( 
u  dx  ^  4  '  ' 


dL 

dr 


_K_ 

pc 


16o1  J 


3k 


i>p( 


dT 

dr 


(34) 


The  universal  velocity  profile  is  expressed  in  terms  of  the  dis¬ 
tance  from  the  tube  wall.  Thus,  a  different  coordinate  system,  winch 
originates  at  the  wall,  must  be  defined.  The  new  coordinate  is 


y  = 


(35) 


Equation  34,  written  in  terms  of  the  new  coordinate,  is 
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u 


dl 

W 


1  d  r  (r  -y) (a+16oT3  +  e)  dT  -i 
Tr'-yT  dy  *•  3kppc  dy  J 


(36) 


The  temperature  gradient  in  the  axial  direction  can  be  determined 
by  performance  of  an  energy  balance  on  a  volume  of  gas  in  the  tube. 

The  temperature  gradient  is  expressed  in  terms  of  the  constant  wall 
heat  flux. 


£L  ,  (37) 

dx  rQpcv 

where  q  >  0  indicates  that  heat  is  transferred  from  the  fluid  to  the 
wall .  0 

Equation  36  may  be  integrated  once  with  respect  to  y.  The  constant 
of  integration  is  determined  from  the  condition  that  the  radial  tempera¬ 
ture  gradient  is  zero  at  the  tube  centerline.  The  results,  when  combined 
with  Equation  37,  are 

r  (rn~yhVl6gT3  +  c),  dX  _  2qo 

L  0  3kRpc  J  dy  “  rQpcv 


y 

u(rQ-y)dy  ]  (38) 


ro 

[  u(rQ-y)dy 

. 

0 


The  dimensionless  variables,  Equations  17,  18,  and  19,  are  now 
introduced  into  Equation  38.  The  dimensionless  form  of  the  energy 
equation  becomes 


*-y*)( 


1  «d03 

Pr  +  PpfRiJ  + 


£ 

V 


V 

u*(rQ*-y*)dy* 


y* 

u*('"0*-y*)dy^  ]  (39) 


where 
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Kays13  gives  the  following  values  for  the  friction  factors  for  fully 
developed  turbulent  flow  in  a  smooth  tube: 

f  =  0.079  Red"0,25  ,  5000  <Red  <  30,000  (41) 


and 

f  =  0.046  Red"0,2  ,  30,000  <  Red  <  1,000,000  (42) 


Either  Equation  41  or  42  may  be  substituted  into  Equation  39,  dependent 
upon  the  Reynolds  number  of  interest. 

The  relations  for  u*  and  e/v  for  the  laminar  sublayer  and  the  buf¬ 
fer  layer  are  again  given  by  Equations  25  and  26.  The  following  values 
were  used  for  the  turbulent  layer,  30  <  y*  <  rQ*: 

u*  =  5.5  +  2.5  In  y*  ,  J  ^  (1  -  ^)  (43) 


The  expressions  for  the  dimensionless  velocity,  u*,  and  e/v  are  substi¬ 
tuted  into  the  integrals  on  the  right  side  of  Equation  39.  When  the 
integrals  are  evaluated,  algebraic  functions  of  rQ*  and  y*  are  obtained. 

The  differential  equations  for  the  layers  are  found  as 
Laminar  Sublayer,  ^  =  F}(e,y*) 

Buffer  Layer,  =  F2(e,y*) 

Turbulent  Layer,  =  F3(e,y*)  (44) 
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where 


Fi(e,y*}  =  {— 


Re/?>"  (ro*-y*>  (p7+~~> 


~}  X 


Pr4Re^ 


{  Ii(r0*)  +  I2(r0*)  +  I3(r0*)-y*2(-y-  -  £■  )  >  (45) 


F2.(0,y*)  =  { 


fed(|)^  (r0*-y*)  (pf  +  7777+  f*  -1) 


}  x 


Pr4ReJ  5 


{  hfr0*)  +  I3(rQ*)  +  3.05  (rQ*y*  -  -  5rQ*  +12.5) 


-  5r0*(y*lny*  -  y*  -  5  In  5  +  5) 


+  S{^~  In  y*  -  -  12.5  In  5  +  6.25))  (46) 


F3(0.y*)  =  { 


Re  H(£f  (r  *-y*)[J-  +  +  J$d  -  j£u 

d  2  0  Pr  Pr4Re3  ^  V 

<  l3(r0*>  "  5,5  (ro*y*  "  "  30ro*  4  45°) 

-  2.5  [r  *(y*lny*  -  y*  -  30  In  30  +  30)] 


}  X 


. v*2  u*2 

f  2.5  (*y-  In  y*  -  *7-  -  450  In  30  +  225)} 


(47). 


and 


MV*)  =  j  U’!(r0*-y*)dy*  =  25(— ^ —  -  |) 


(48) 


! 


I 


30 


I2(r0*) 


u*(r0*-y*,)dy*  =  ••  3.05  (25rQ*,  -  437.5); 


+  5ro*(30  In  30  -5  In, 5  -  25) 


-.5(450  In  30  -  12.5  In  .5  218.75) 


(*9) 


o  ! 


*2 


1 3 ( r0*)  ='  u*.(rQ*-y*)dy*  =  5.‘ 5  (-| - ,30ro*  +  450) 


30 


+  2.5  [r  *(jr  *lnr  *  -  30  In  30  -  rQ*  +  30)] 


s  r 


*2 


*2 


-  2.5  Hj—  In  rQ*  -  -  450  In  ,30  +  225] 


(50) 


*  i  * 

The  complete  temperature  distribution  in  the  tube  can  be  determined  from 
the  solution  of  these  equations. 

1  The  value  of  r  *,  which  is. a  function  of  the  Reynolds  number,  must 
be  known  before  theaifferential  equations  can  be  solved.  The  mass  flow 
rate  in  the  tube  >i s  >  i 


! 

’t 


m  =  2tip  j  urdr 
o 


(51) 


The  Reynolds  number  may ; be  defined  in  terms  of  the  ma,ss  flow  rate  as 


Re 


2m 


d  ,  npvr. 


(52) 


When  Equations  51  and  52  are  combined  arid  expressed  in  dimensionless 


form,  the  following  expression  may  be  obtained  for  rQ*: 


!  *  4 

r  *  = 

o  Re. 


V 


u*(rQ*-Y*),dy* 


(53) 
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or 


*.  A 


EWV)  +  I2<’V>  +  M^*)] 

Cl 


Equation  54  may  be  solved  numerically  for  r  *  by  the  Newton-Raphson 
method.14 

The  Mussel t  number  can  be  calculated  when  the  complete  temperature 
distribution  is  known.  The  Mussel t  number  is  defined  as 

...  _  hd  _  V 

Nu "  t  "  K[r^r“r  (55) 

'  CO  o 


The  difference  between  the  wall  temperature  and  the  centerline  tempera¬ 
ture,  or  the  difference  between  the  wall  and  bulk  temperature,  may  be 
used  to  evaluate  Equation  55.  Introduction  of  the  dimensionless  variables 
into  Equation  55  gives  the  Nusselt  number  on  the  basis  of  maximum  tempera¬ 
ture  difference. 

Re.  Pr 

Nuc  =  IF^T  (55) 


which  can  be  evaluated  when  the  dimensionless  centerline  temperature  has 
been  calculated. 

The  bulk  temperature  is  defined  as 


'b“CT 


uTrdr 


The  dimensionless  variables  are  introduced  into  Equation  57  which  is 
then  substituted  into  the  dimensionless  form  of  Equation  55.  The  fol¬ 
lowing  expression  is  obtained  for  the  Mussel t  number  on  the  basis  of 
the  bulk  temperature  difference. 


Nub  = 


Rej  Pr  r  *2 
d  o 


» 

ffi  u*e(r  *-y*)dy*  -  r  *2e. 

0  Or 
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Equation  58  v.ars  be  numerically  evaluated  when  the  temperature  disttibu- 
tion  is  known. 

The  governing  equations  were  evaluated  for  the  case  where  Cj  is 
specified,  and  -for  the  case  where  the  wall  and  centerline  temperature 
were  specified.  An  iterative  solution  was  required  in  the  latter  case. 
The  analytical  results  are  presented  in  the  following  section. 

ANALYTICAL  RESULTS 


The  differential  equations,  derived  in  the  previous  section,  were 
numerically  evaluated  by  use  of  the  fourth  order  Runpe-Kutta  algorithm. 
The  results  for  turbulent  flow  over  a  flat  surface  are  shown  in  Figures 
2,  3,  and  4.  The  Nusselt  number  :s  plotted  against  the  Reynolds  number 
as  a  function  of  Sx  and  the  dimensionless  wall  temperature  in  Figure  2. 
Radiation  effects  are  important  if  Re,{  is  less  than  2  x  105,  For  laroer 
values  of  the  Reynolds  number,  turbulent  convection  is  the  dominant 
mode  of  heat  transfer.  The  relative  importance  of  radiation  heat  trans¬ 
fer  is  strongly  dependent  on  the  value  of  Sx*  Radiation  is  important 
when  5X  is  large.  The  value  of  £x  is  large  when  the  temperature  is 
high  or  the  absorption  coefficient  is  small.  The  value  of  the  radiation 
contribution  also  depends  on  the  value  of  the  bore  surface  temperature. 

A  higher  surface  temperature  results  in  a  relatively  laraer  heat  trans¬ 
fer  contribution  by  radiation.  The  Nusselt  number  is  plotted  against 
the  Reynolds  number,  in  Figures  3  and  4,  as  a  function  of  centerline 
temperature  and  the  parameter  u  /lev .  A  bore  surface  temperature  of 
1000°R  is  shown  in  Figure  3,  and  a  bore  surface  temperature  of  2000°R 
is  shown  in  Figure  4.  The  radiation  contribution  is  important  when 
Re  is  less  than  2  x  105.  Radiation  is  significant  at  the  higher 
Reynolds  numbers  only  if  the  parameter  uro/kv  is  large.  This  again 
means  that  the  absorption  coefficient  must  be  small.  The  values  of 
uOT  and  v  are  known  from  the  flow  conditions  in  the  gun  tubes.  Thus, 
the  relative  importance  of  radiation  can  be  quickly  estimated  when  the 
absorption  coefficient  is  known.  The  radiation  contribution  for  a 
given  centerline  temperature  is  greater  when  the  surface  temperature 
is  higher  as  shown  in  Figures  3  and  4.  This  means  that,  as  the 
barrel  temperature  rises,  the  heat  flux  does  not  decrease  as  much  as 
would  be  calculated  by  the  convective  heat  transfer  theory. 

The  Nusselt  number  for  fully  developed  turbulent  flow  is  plotted 
in  Figures  5,  6,  and  7.  The  Nusselt  number  is  plotted  against  the 
Reynolds  number  as  a  function  of  ?x  and  the  dimension less  wall  tempera¬ 
ture.  Radiation  is  important  when  the  Reynolds  number  is  less  than 
2  x  105.  Turbulent  convection  heat  transfer  is  dominant  for  larger 
Reynolds  numbers.  The  relative  importance  of  radiation  is  stronnly 
dependent  upon  the  value  of  Radiation  is  important  for  large 
values  of  £<]>  which  implies  that  the  temperature  must  be  hi  oh  or  the 
absorption  coefficient  small.  Higher  bore  surface  temperatures  result 
in  a  relatively  higher  contribution  of  radiation  to  the  total  heat 
transfer.  The  Nusselt  number  is  plotted,  against  the  Reynolds  number 
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FIGURE  4.  NUSSELT  NUMBER  VERSUS  REYNOLDS’  NUMBER -AS  A  FUNCTION  OF 

/kv  FOR  TURBULENT  FLOW  OVER  A  FLAT  PLATE  WITH 
THICK  RADIATION 

!  •  I 


Tg  AND  u„ 
OPTICALLY 
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FIGURE  6.  NUSSELT  NUMBER,  BASED  ON  CENTERLINE  TEMPERATURE  DIFFERENCE, 
VERSUS  REYNOLDS  NUMBER  AS  A  FUNCTION  OF  OPTICAL  DIAMETER 
FOR  OPTICALLY  THICK,  FULLY  DEVELOPED,  TURBULENT  FLOW  IN  A 
TUBE.  Tg  =  6000°R,  T0  =  1000°R 
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FIGURE  7.  NUSSELT  NUMBER,  BASED  ON  CENTERLINE  TEMPERATURE  DIFFERENCE, 
VERSUS  REYNOLDS  NUMBER,  AS  A  FUNCTION  OF  OPTICAL  DIAMETER 
FOR  OPTICALLY  THICK,  FULLY  DEVELOPED  TURBULENT  FLOW  IN  A 
TUBE.  T a  =  6000°R,  Tn  =  2000°R 


as  a  function  of  the  optical  diameter,  as  shown  in  Figures  6  and  7.  The 
centerline  temperature  is  6000°Rin  both  of  these  figures:  the  surface 
temperature  is  1000°R  in  Figure  6,  and  2000°R  in  Figure  7.  The  contri¬ 
bution  of  radiation  is  a  strong  function  of  the  optical  diameter.  For 
large  values  of  the  optical  diameter  (Td  >  1000),  the  effects  of  radia¬ 
tion  are  small.  This  occurs  because  the  major  portion  of  the  radiation 
emitted  by  the  gas  is  absorbed  before  the  radiation  strikes  the  bore 
surface.  The  effects  of  radiation  are  greater  for  higher  wall  tempera¬ 
tures. 

The  data  presented  in  the  figures  may  be  used  to  quickly  estimate 
the  importance  of  radiation  heat  transfer  in  any  gun  tube  in  which  the 
flow  conditions  are  known.  The  XM140  gun  tube  may  be  considered  as  an 
example.  The  projectile  velocity,  gas  density,  pressure,  and  gas  tem¬ 
perature  have  been  plotted  against  projectile  displacement  by  Dahm  and 
Anderson.10  When  the  base  of  the  projectile  is  three  inches  from  the 
breech  end,  the  conditions  in  the  gun  tube  are 


L  =  3  in 

P  =  25,000  lb/in2 

Vp  =  325  ft/sec 

p  =  7  lb/ftj 

T  =  5000°  f< 


The  value  of  the  viscosity  can  be  calculated  from  the  relation10 


T  o.74 

"  *  °-044  <S3o> 


lb 

hr  ft 


(59) 


The  value  of  the  Rossela"d  mean  absorption  coefficient  can  be  estimated 
from  Equation  10.  These  data  were  used  to  calculate  Rex  and  u^/k^v  as 

functions  of  x.  The  calculated  values  are  listed  below: 


XJjn). 

!fx 

0 

0 

0.5 

2.34  x 

1.0 

9.36  x 

1.5 

2.11  x 

2.0 

3.75  x 

2.5 

5.85  x 

3.0 

8.43  x 

0 

10 

0.55  x 

102 

105 

1.10  x 

102 

106 

1 .65  x 

10* 

106 

2.20  x 

102 

106 

2.75  x 

102 

!0b 

3.30  x 

102 
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The  boundary  layer  thickness  can  be  calculated  from13 


6  =  0.3814  x  Re  "°-2 

A 


(60) 


The  maximum  boundary  layer  thickness  is  0.047  inch  which  is  much  less 
than  that  of  the  tube  radius,  0.59  inch.  The  flow  is  not  fully  devel¬ 
oped,  and  the  effects  of  radiation  can  be  estimated  by  reference  to 

Figures  3  and  4.  Radiation  is  an  important  mode  of  heat  transfer  in 

the  first  two  inches  of  the  gun  tube  where  the  Reynolds  number  is  less 

than  106.  At  higher  values  of  the  Reynolds  number,  convection  is  pre¬ 

dominant.  The  effects  of  radiation  are  small  because  of  the  large 
value  of  the  absorption  coefficient  at  this  pressure.  Nearly  all  of  the 
energy  emitted  by  a  volume  of  gas  is  absorbed  by  other  gas  volumes  before 
it  reaches  the  bore  surface. 

When  the  base  of  the  projectile  is  12.5  inches  from  the  breech  end, 
the  conditions  in  the  gun  tube  are 


L  =  12.5  in 

P  =  12000  psi 

Vp  =  1750  ft/sec 

p  =  4.5  lb/ft3 

T  =  4000°R 


These  data  were  used  to  compute  the  values  of  Reu  and  u  /kDv  qiven  below. 

X  m  K 


;  (in) 

Rex 

£ 

0 

0 

0 

1.0 

9.61 

X 

105 

2.34 

X 

I02 

2.0 

3.84 

X 

106 

4.70 

X 

102 

3.0 

8.64 

X 

108 

7.06 

X 

102 

4.0 

1.54 

X 

107 

9.45 

X 

102 

5.0 

2.40 

X 

107 

1.17 

X 

103 

6.0 

3.46 

X 

1°7 

1.41 

X 

103 

7.0 

4.71 

X 

107 

1.65 

X 

103 

8.0 

6.15 

X 

107 

1.88 

X 

103 

9.0 

7.79 

X 

107 

2.12 

X 

103 

10.0 

9.61 

X 

107 

2.36 

X 

103 

11.0 

1.17 

X 

108 

2.58 

X 

103 

12.0 

1.38 

X 

108 

2.82 

X 

103 

12.5 

1.50 

X 

10s 

2.94 

X 

103 
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The  boundary  layer  thickness  calculated  from  Equation  60  is  0  009  inch 
which  is  much  less  than  that  of  the  tube  radius.  The  effects  of  radia¬ 
tion  can  again  be  estimated  from  Figures  3  and  4.  Radiation  is  signifi¬ 
cant  in  the  first  three  inches  from  the  breech.  Reycnd  this  distance, 
turbulent  convection  is  dominant. 

In  both  of  the  examples  given  above,  radiation  was  important  only 
in  the  region  near  the  breech.  This  will  be  true  for  ail  small  caliber 
weapons  because  the  propellant  gas  pressures  are  very  high.  At  these 
high  pressures,  the  gas  is  very  dense  and  is  nearly  opaque  to  thermal 
radiation.  Thus,  while  a  large  quantity  of  radiant  energy  is  emitted 
by  the  propellant  gas,  most  of  this  energy  is  absorbed  before  it  strikes 
the  bore  surface.  If  the  projectile  velocity  and  the  tube  pressure  were 
lower,  radiation  would  be  a  more  important  mode  of  heat  transfer.  TTie 
absorption  coefficient  would  be  smaller,  and  more  radiation  would  -'trike 
the  bore  surface.  Also,  the  turbulent  convective  heat  transfer  would  be 
less.  Thus,  radiation  heat  transfer  would  be  relatively  more  important. 
These  conditions  may  exist  in  some  artillery  weapons. 

SUMMARY 

Thermal  radiation  in  gun  tubes  was  analytically  studied.  A  simpli¬ 
fied  physical  model  was  used  in  this  analysis  With  this  mode!,  steady 
turbulent  flow  of  an  optically  thick  gray  gas  in  the  gun  tube  war.  con¬ 
sidered.  The  flow  in  both  undeveloped  and  fully  developed  turbulence  was 
considered.  The  governing  equations  were  derived  and  proa rammed  for 
numerical  evaluation  The  numerical  results  were  presented  m  dimension¬ 
less  form  so  that  they  could  be  used  to  estimate  the  contribution  oi 
radiative  heat  transfer  in  any  gun  tube  for  which  the  iiow  v-onditionb 
are  known.  The  contribution  or  the>mai  radiation  heat  transfer  in  the 
XM140  gun  tube  was  investigated,  as  cm  example.  Radiation  was  ’mpo  ■ 
tant  near  the  breech  end.  At  locations  cr  more  tnan  a  few  mules  worn 
the  breech  end,  turbulent  convection  was  tne  dominant  mode  of  heat 
transfer.  The  primary  reason  that  radiation  is  insignificant  at  dis¬ 
tances  further  from  the  breech  is  that  the  absorption  coefficient  oi  the 
propellant  gas  is  very  large  Whiie  a  significant  amount  of  rad; am 
energy  is  emitted  by  the  gas,  most  of  the  radiant  energy  is  aoso.-oed 
before  it  strikes  the  bore  surface.  Radiation  would  De  mo*e  mpj'tant 
in  guns  operating  at  lower  pressures  and  having  tower  p- eject’ To  veloc¬ 
ities. 
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